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Abstract 

We employ a three fluid model in order to construct a cosmological model in the Friedmann Robertson 
Walker flat spacetime, which contains three types of matter dark energy, dark matter and a perfect fluid 
with a linear equation of state. Dark matter is described by dust and dark energy with a scalar field with 
potential V{<j)). In order to fix the scalar field potential we demand Lie symmetry invariance of the field 
equations, which is a model-independent assumption. The requirement of an extra Lie symmetry selects 
the exponential scalar field potential. The further requirement that the analytic solution is invariant under 
the point transformation generated by the Lie symmetry eliminates dark matter and leads to a quintessence 
and a phantom cosmological model containing a perfect fluid and a scalar field. Next we assume that 
the Lagrangian of the system admits an extra Noether symmetry. This new assumption selects the scalar 
field potential to be exponential and forces the perfect fiuid to be stiff. Furthermore the existence of the 
Noether integral allows for the integration of the dynamical equations. We find new analytic solutions to 
quintessence and phantom cosmologies which contain all three fiuids. Using these solutions one is able to 
compute analytically all main cosmological functions, such as the scale factor, the scalar field, the Hubble 
expansion rate, the deceleration parameter etc. 

PACS numbers: 98.80.-k, 95.35.+d, 95.36.+X 
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1 Introduction 

The recent cosmological data indicate that the universe (a) is spatiaUy flat, (b) has suffered two acceleration 
phases. An early acceleration phase (inflation), which occurred prior to the radiation dominated era and a 
recently initiated accelerated expansion. The source for the late time cosmic acceleration has been attributed 
to an unidentified type of matter, the dark energy (DE). DE contrary to the ordinary baryonic matter has a 
negative pressure which counteracts the gravitational force and leads to the observed accelerated expansion. 
The nature of DE is still an open question. 

The simplest DE probe is the cosmological constant A (vacuum) leading to the ACDM cosmology [3 [2J 13] • 
However, it has been shown that ACDM cosmology suffers from two major drawbacks known as the fine 
tuning problem and the coincidence problem [4]. Besides ACDM cosmology, many other candidates have 
been proposed in the literature, such as time-varying h.(t) cosmologies, quintessence, fc— essence, tachyons, 
modifications of gravity, Chaplygin gas and others [SJ IH [71 HI HI ■ 

In addition to DE and the ordinary baryonic matter, it is believed that the Universe contains a third type of 
matter, the dark matter (DM). This type of matter is assumed to be pressureless (non-relativistic) and interact 
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very weakly with the standard baryonic matter. Therefore its presence is mainly inferred form gravitational 
effects on visible matter. 

In the following we consider a model of the Universe which contains three types of matter. The DM is 
modeled by a dust fluid, the DE by a scalar field and the rest of matter by a perfect fluid. Since a scalar field 
can be considered to be a perfect fluid (see below) in effect we model the matter in the Universe in terms of two 
perfect fluids and a dust. All three fluids are assumed to be self-interacting and minimally coupled to gravity. 

The problem with the above scenario is that there does not exist an underline principle which will specify 
uniquely the potential 1^(0) of the scalar field. Indeed in the literature one finds several potentials V{(j)) such 
as exponential, power law, hyperbolic etc. Therefore it is not possible that one could find an analytic solution 
of this model even if an EoS has been assumed. 

The aim of the present work is twofold (a) to propose a geometric principle ('selection rule') for specifying 
the potential V{(j)) and (b) to solve analytically the system of the resulting field equations. Concerning (a), 
following a recent paper [TT], we propose that the potential should be selected by the requirement that the 
dynamical system of the three fluids admits an additional Lie or Noether symmetry. This point of view has also 
been considered in [H [H [H [H] . 

The main reason for the consideration of this hypothesis is that the Lie/Noether point symmetries provide 
first/ Noether integrals, which assist the integrability of the system. A fundamental approach to derive the 
Lie point and the Noether symmetries of a given dynamical system moving in a Riemannian space has been 
proposed recently in [16] . A similar analysis can be found in [171 Ull Ull IlHl [21] . 

Concerning the analytic solutions we restrict our considerations to a flat FRW spacetime. Although there 
is a great number of papers devoted to the dynamics of scalar field minimally coupled to matter, not much is 
known about analytic exact solutions of these models. Most of the solutions correspond to spatially fiat FRW 
models with no other source but the scalar field [22l [23l [24l [25l [26] . Recently in [11] the analytic solution for 
matter in the form of dust and a scalar field has been given. A solution with two scalar fields (where one of 
them is a kination i.e. stiff matter) and an exponential potential is given in [27]. 

Few exact solutions are known with spatial curvature [28j [29] . Even less solutions are known for a perfect 
fluid and a scalar field [SJ [301 [SIl [32] ■ In particular in [35] the authors consider a prefect fiuid minimally 
interacting with a scalar field having an arbitrary potential V{(j)) and, by considering as variable the scale factor 
a, they develop an approach which replaces the potential V{(j)) with an auxiliary arbitrary function (F(a) in 
their notation) and a constant (C) and give the implicit solution of the problem. They consider an arbitrary 
specification of the function F(a) and the constant C, which leads to an exponential potential and by (making 
some further assumptions) they find a set of analytic solutions in terms of the adiabatic index 7 of the fluid. 
In a different approach to the same problem in |12| . the authors employ ad hoc an exponential potential and 
assume that the resulting system of equations admits a Noether symmetry. They are not able to find an analytic 
solution and continue with numerical solutions. This is not necessary according to the results of [32] and our 
results below. 

In [331 134] the authors in an attempt to unify the early time (phantom) inflation with the late time (phantom 
or not) acceleration derive specific solutions within which both the early and the late phases of evolution of the 
universe occur. The proposed scenario is based in the introduction of an additional unspecified function io{(j)) 
in the kinetic term of the scalar field Lagrangian. Although this function can be absorbed into the field during 
inflation and acceleration by a redeflnition of the scalar field, this is not possible right on the transition point, 
at which this term connects the two phases of evolution smoothly. A similar scenario considers two scalar fields. 

There does not appear to exist a solution concerning a perfect fiuid, a dust and a scalar field all minimally 
coupled to gravity and non-interacting in a spatially fiat FRW spacetime. In the following using symmetry 
assumptions only, we obtain many of the above solutions and also a new solution which involves all three fluids. 
This solution contains the solutions of [27] as a special case. 

The structure of the paper is as follows. The basic equations and the resulting dynamical system are 
presented in section 2. In section 3 we review briefly the basics of the theory of Lie and Noether symmetries. 
In section 4 we define the scalar field potential and subsequently we determine analytic solutions of the field 
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equations using an extra Lie symmetry and the Lie invariance of the solutions with respect to this symmetry. 
In section 5 we determine a new solution of the three fluid dynamical system assuming only an extra Noether 
symmetry. The main conclusions are summarized in section 6. 



2 The dynamical system 

We consider the flat FRW spacetime in canonical coordinates: 

ds^ = -dt^ + (t) {dx^ + dy^ + dz^) (1) 

with comoving observers u"" = S^. We also consider that spacetime contains dark energy, dark matter and 'other 
matter'. DM is assumed to be described by dust {pdm = 0) with energy momentum tensor 

DAlTab = flDMUaUb- (2) 

The 'other matter' is assumed to be described by a perfect fluid with energy momentum tensor 

PpTab = ^^B'^aUb + psKb (3) 

and a linear EoS 

PB = h-l)^iB , i<7<2 (4) 

where 7 is the adiabatic index of the fluid and hab ~ gab+UaUb is the tensor projecting normal to the four- velocity 
It". For 7=1 the fluid is dust; for 7 = 2 is stiff matter and for 7 = | is a radiation fluid. 

Finally DE is described by a scalar field 4> (quintessence or phantom) rolling down a potential V{4>) whose 
energy momentum tensor is 

c^Tab = Qe^' + V{(l))^ UaUb + Qef/)' - hab. (5) 

When £ = +1 we have quintessence (real (p) and when e = — 1 we have a phantom field (imaginary (p). 
It follows that (for the comoving observers) the scalar field may be considered as a perfect fluid with energy 
density and isotropic pressure: 

H = ^e<^' + ^(<^) (6) 

P^ = \e^'' -V{^). (7) 

Equations dH),© imply the EoS: 

^Ji^-p^ = 2V{cj>). (8) 

from which follows that the definition/specification of a potential function for the scalar fleld (f){t) essentially 
defines an EoS for the perfect fiuid defined by the scalar field. It also follows that stiff matter may be considered 
as a homogeneous, massless, free (i.e. with vanishing potential) scalar field minimally coupled to gravity. Because 
■2^ = 1 this scalar field is a kination. 

All three fluids are assumed to interact minimally, which implies the three independent conservation equa- 
tions: 

z5A/n' = 0, PFT'^b^^O, 0Tt = O. (9) 

There are six unknowns in the problem i.e. the a(t),/i^, pB, fJ-uMT^'it): ^('/') ^-nd five field equations i.e. 
Einstein's equation: 

Gab —<t> Tab +DM Tab +PF Tab (10) 

the three conservation equations ([9|) (equation 4,T°:l = is the Klein Gordon equation) and the EoS (jlj. We 
need one more equation which must be given by an independent assumption. 
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However there is not a general consensus or accepted practice as to the nature of this assumption. In the 
literature one finds the following types of assumptions/equations: 

(i) Equations which are based on the dynamical system approach. That is, in the case the field equations follow 
from a Lagrangian one demands that the Lagrangian admits a Noether symmetry. This assumption results in an 
extra Noether integral, which provides the required extra equation. With this approach there have been given 
the analytic solutions for various probes (UCDM, exponential) [TTJ [T^l [H] , and in other models e.g. modified 
gravity [351 ISSl [3Z1 [SE [Ml US HIl HI] • An alternative but similar approach is to demand that the field equations 
admit an extra Lie (not necessarily a Noether) symmetry and require that the solution is invariant under the 
action of the Lie symmetry, (ii) Equations which are defined by an ansatz on the scale factor [24l [25l [26] . (iii) 
Equations which are ad hoc statements on the scalar field potential. 

In the present paper we follow the dynamical system approach and demand that the system admits an 
extra Lie symmetry and subsequently an extra Noether symmetry. The Lie symmetry requirement specifies the 
potential to be exponential and produces (with a additional invariance assumption) the perfect fluid solutions 
found in [32| and the empty space solution of [24l [26] . The requirement of Noether symmetry in addition of 
fixing the potential, it constraints the perfect fluid to be stiff and determines the analytic solution of the three 
fluid problem. 

2.1 The field equations 

Einstein field equations for the metric ([1]), comoving observers and a perfect fiuid of energy density /i^ and 
isotropic pressure ps minimally coupled to a scalar field (p with scalar field potential V{(l)) and DM (dust) with 
energy density are: 

k 

= ^ O^DM + f^B + M0) (Hubble equation) (11) 
a 1 k 

— h -H^ = — — {pb +Pch) (deceleration equation) (12) 
a 2 2 

where k = SttG is Einstein's gravitational constant. The deceleration equation (jl2p can be written in the 
alternative forms 

H+^H^^-^ipB+p^) (13) 



H 



^ [t^B + PB + e4>') ^ 0. (14) 
Equations (|lip . (IT2l) are supplemented by the Klein Gordon equation 



(l,^~3H(j)-eV^ (15) 

and the conservation equations: 

AB + 3(/iB+Pij)if = (16) 

^DB + ''^P'DbH = 0. (17) 

The conservation equation ([16)) and the EoS ([4]) imply 

^lB = C^a"'^ , PS = (7 - 1) , K 7 < 2 (18) 

where C ^ is an integration constant. For DM the conservation equation (IT71) gives: 

E , . 

^^DM = — 19) 

where E is an integration constant. Replacing these results in the field equations pT|) . p2)) , (ITSt we obtain the 
system of equations 

Saa^ - l-eka^^^ - ka^V (0) - kCa'^'^''-^^ = E (20) 
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a+—a^ + -ea^^ - -aV + - {-f ~ 1) Ca'-''^ 
2a A 2 2 

3 • 

(j> + -d(j) + eVA, 
a 

where E = kE. Equation (j20p is written in the equivalent form 

i/2 = I (^Ea-^ + + ie^' + ^(0)^ . (23) 

2.2 Looking for analytic solutions 

In order to find an analytic solution of the system of equations (pil)) - (|^^ we need either two independent 
assumptions, which will result in two independent conditions, or one stronger assumption which will result in 
two conditions. In the following wc formulate these assumptions in terms of Lie and Noether symmetries. 

Assuming the unknown potential 1^(0) to be an externally specified function, the differential equations (j2ip . 
(1221) contain only the unknowns a{i), (/'(t). Hence it is possible to consider a new autonomous two dimensional 
dynamical system with variables {a{t),4'{t)} defined by equations (PT|) . ([^ and treat the third field equation 
(l20l) (equivalently (1231) ) as a constraint. Following this remark we look for a Lagrangian producing equations 
(f2T|) . ([22|) . It is easily seen that this Lagrangian is 

L = (^aa^ ~ hiea^'cjf^ + ka^V (0) + kCa'^^^-'^'i (24) 

with Hamiltonian 

E = 3ad2 - -fcea^0^ - ka^V (0) - fcCa'^^'^-^). (25) 

We note that the Hamiltonian is the constraint condition (|20|) . that is i? = const. It is possible to give a 
geometric interpretation of the Hamiltonian E. Indeed the dynamical system defined by the Lagrangian ([M]) is 
autonomous hence admits the Noether symmetry dt whose Noether integral is precisely the Hamiltonian. We 
note that E — only when DM is absent. We arrive at the conclusion: 

The dynamical system defined by the three fluids considered above is equivalent to a new autonomous two 
dimensional dynamical system {a{t) , (l){t)} with Lagrangian {24-^ whose Hamiltonian E is a constant which 
vanishes iff ^ij^j^j = 0. 

In the following we use the Lie and the Noether symmetries of the two dimensional dynamical system in 
order to produce analytic solutions. 



(21) 
(22) 



3 Lie and Noether symmetries 

Before we proceed we review briefly the basic definitions concerning Lie and Noether symmetries of systems of 
second order ordinary differential equations (ODEs) 

^ijj' {t,x^ ,x') . (26) 

A vector field X ~ ^{t, x^) dt + rf {t, x) di in the augmented space {t, a;*} is the generator of a Lie point 
symmetry of the system of ODEs ([26|l if the following condition is satisfied [ITl [18] 

A:[21 (i*-u;(t,a;^i^)) =0 (27) 

where ATI^I is the second prolongation of X defined by the formula 

= S^dt + rfd, + [ij' - i'^) d^. + {if - x'l - 2i'e) d^. . (28) 
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Condition (|27p is equivalent to the relation: 



= A(a;°) A (29) 

where X'^' is the first prolongation of X and A is the Haniiltonian vector field: 

A = dt + xd^+ uj' (t, , x^) di, . (30) 

If the system of ODEs results from a first order Lagrangian L = L (t, x^ , i-') , then a Lie symmetry X of the 
system is a Noether symmetry of the Lagrangian if the additional condition is satisfied 

XWL + L^^f (31) 
dt dt ^ ' 

where / = f (t,x^) is the gauge function. To every Noether symmetry there corresponds a first integral (a 
Noether integral) of the system of equations I^E^ which is given by the formula: 

where Eh is the Hamiltonian of L. 

The vector field X for the Lagrangian ([24|) is 



X = £,{t, a, (j)) dt + rj^ {t, a, cj)) da + rj^ {t, a, (f) (33) 

and the first prolongation 

^''1 = + Vada + V^d^ + {va ^ a^) da + (?70 - 0e) d^. (34) 

Having given the basic formulae for the Lie and Noether symmetries we look for analytic solutions of the 
system of equations (PO]) - (1^ using Lie and Noether symmetries. 

To simplify the calculations we introduce a new variable r with the relation a — r's . In the new variables 
r, (p the Lagrangian (|24|) takes the form 



i = i _ skr^4f^ + kr^V (</>) + kCr-^^^-^\ 



(35) 



For a general V{4)) this Lagrangian admits only the standard Lie/Noether symmetry dt- We are looking 
for the possibility of extra Lie and Noether symmetries for special forms of the potential V{4)). To find these 
potentials we apply the results of [16j , which concern all autonomous two dimensional dynamical systems moving 
in flat (Euclidian or Minkowskian) space and admit Lie and/or Noether symmetries. 

In order to apply the results of |16j we need a fiat metric and a potential function for the dynamical system. 
This is done by decomposing the Lagrangian psp into kinetic energy, which defines the kinetic energy metric, 
and potential energy which defines the potential. We consider the metric to be: 

dsl = - ekr^ip^ (36) 

and the potential function: 

H/(r,0) = -fcrV ((/))- /fcCr-2(7-i). (37) 

It can be shown that the Ricci scalar of the metric ds\ vanishes, hence the space {r,0} is the (Minkowski 
two dimensional space in polar coordinates). Therefore the results of [16j apply and we make use of the tables 
in that paper to read the appropriate potentials. 
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4 Lie symmetry 



We start with the demand that the system of equations ([?T|) . admits an extra Lie symmetry (apart of 
dt). From Table 10 Line 5 of [IS] we find that the system of equations admits an extra Lie symmetry for the 
potential V {(p) = Voe~'^'^ where Vq , d are constants with d 0; Furthermore the Lie symmetry vector is 



or, in the original coordinates a{t), cj) (t): 



X = jtdt + rdr + (38) 



For the scalar field potential V{(j)) = Vqc '^'^ the field equations become: 

(^3ad^ - ifcea^^^^ - ka^V {(f>) ~ kCa-^^^-^^ = E (40) 
a + + -^eacp^ - ^aV^e"'^ + Ca'-^'' = (41) 



'<j) + 3-^ - deVoe-'^'l' = 0. (42) 

where 

Mb = Ca-^'' , PB^il-l) Ca-^-*. (43) 

= l^i'' + y^^-"^ , = 1^4'' - Voe-'t (44) 
We look for analytic solutions of these equations. 

4.1 Solution using Lie invariance 

To find a solution of the system of equations we need one further condition. We require the solution to be 
invariant under the action of the extra Lie symmetry pop . To find this solution we compute the zeroth order 
invariants of X using the associated Lagrange system (recall that 1 < 7 < 2 and d ^ 0): 



dt da 

l« IT 

The solution of the system is the two parameter family of functions: 



7t ia ^ ■ 



2 2 

a(t) = t— , (j){t) ^ - Int. (45) 
d 

Subsequently we demand that the (invariant) a{t), (j){t) we computed are solutions of the field equations 
(PTjl . (HH). This fixes the constants Vo, C in terms of 7, d as follows: 

_ 2£(2-7) 4(d^-3£fc7) 

We note that in both solutions the DM is eliminated because i? = 0. Concerning the properties of these 
solutions we have the following. 

A. In the quintessence case the scalar field ((> is real. Because 7 < 2, Vq > the V{4)) > 0, therefore the 
potential is "repulsive". The constant C is related to the matter energy /i^ via relation P5|) . For the perfect 
fluid fj,g and C are positive which requires that d^ > 3k^. 

B. In the case of phantom field e — —1, (j) is complex, d is complex but V {(p) is real. Furthermore Vq > 0, 
and C ^ provided Id^l ^ 3^7. 
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When l^^l = 3fc7 then C — and we have the empty space solution with either a quintessence (e = +1) or 
a phantom field (e = —1). This solution has only the unspecified parameter 7. We note that when 7 = 2 (stiff 
matter) the constant Vq = which is not acceptable. Therefore the case 7 = 2 must be treated separately. 

Perfect fluid solutions with a linear EoS and a scalar field with exponential potential have been considered 
previously in [32] ■ It can be shown that our solution contains the solution of [32 . It is to be noted that we have 
derived the solution using only fundamental symmetry assumptions and not ad hoc statements. Furthermore 
we have covered the case of the phantom field. 

For all cases with 7 < 2 we have a genuine perfect fiuid solution minimally interacting with a scalar field . 
At late time both fi^jPs — >■ 0, that is, the solution p5t for late time tends to the empty space solution with a 
scalar field. This results agree with known results (see [22 l [24 j [26]). 



5 Noether symmetry 

In an alternative approach we continue a stronger symmetry assumption and require that the dynamical system 
defined by the Lagrangian psp admits a Noether symmetry. It is known that, contrary to the Lie symmetry, the 
Noether symmetry gives two conditions on the system of ODEs which suffice to produce the analytic solution 
of the dynamical system 

Because the Noether symmetries are Lie symmetries and we have only one extra Lie symmetry, the Noether 
symmetry must coincide with the Lie symmetry we found above. This leads to the value of the parame- 
ter 7 = 2. Therefore the requirement of Noether symmetry has the following implications: (a) Selects the EoS 
Pb — Pb; (b) the original field equations ((2T]) . (|22]) are non- homogeneous, (c) Selects the scalar field potential 
V{(j)) = V^e^'^'^ where Vq is a constant, (d) fixes the Noether symmetry vector ([55]) and (e) Provides the Noether 
integral: 

4 

I2 = 2tE - Aa^a + -eka^ij). (47) 
a 

In the original coordinates a, (j) we have 
Lagrangian: 

L = (^aa^ - htea!^^^^ + kVoa^e-'^'^' + kCa'^ (48) 

Hamiltonian: 

E = ^3ad2 - hisa^ijf^ - fca^Foe"'*'^ - kCa-^ (49) 

Equations of motion: 

a + ^ ^eadp^ - ^aV^e-"^ + ^Cfl-^ = (50) 
2a 4 2 2 

^ + -a^ - deVoe-'^'f' ^ 0. (51) 
a 

We conclude that the requirement of the existence of an extra Noether symmetry results in a model consisting 
of a kination, a scalar field with an exponential potential and dust. 



5.1 The analytic solution 



We determine the analytic solution of the system of field equations (150]) . ([51]) under the constraint ([49]) . In the 
solution we follow a procedure initiated in [22]. We introduce new variables {u, v} with the relations 



— 0+-ln(a3) 

.^0+lln(a3) 
4^2 ^ ' 



(52) 
(53) 
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In the new variables the Lagrangian (|^5)) is 



where ii' 



■7=, Vb = A:Vb. Next we consider a change in the time coordinate 



dt~\ 4 

Using the variation integral the Lagrangian becomes 

L = e(i-^)"e(i+^)- (u'v' + 1 + C'e-2(i-^)"e-2(i+^)- 

where C* — j^y^C and the Hamiltonian 

E = e(i~^)"e(i+^)'' (i,V - 1 - Ce-2(i-^)«e-2(i+^)- 

The resulting equations of motion are: 

u" + il-K) u'^ - (1 + X) (1 - C'e-2(l-^)«e-2(l+/^)^.^ 
z;" + {1 + K) v'^ -{l-K)(l- (7e-2(i-^)«e-2(i+^)-) 

We consider the cases K — 1 and K ^ 1. 

5.1.1 The case K = 1 

For K = 1 the system of equations ([58| . ((59l) becomes: 



u"-2(l-Ce 



and the Hamiltonian (l5 



The solution of the system §0^ - §^ is 



E^e^" (u'v' -l-Ce- 



u{t) = + -C\n{T + c) + {E + 2c) t 
w (t) = i In (2r + 2c) . 
Replacing in ([5^ . (|53p and ((M)) we find the analytic solution 

where 



(r) = V'2(r + cf (T + c)e-' + (^+2-> 



0(r) 



/6fce 



1 , /I 



+ (i; + 2c) T + - In - (r + c) 



nC-1 



5.1.2 The case K 1 

Subcase K < 1 In this case the analytic solution of the system of equations (ISOMQp is (for details see 
Appendix) 



(r) - - (1 - Kf-"^ (1 + Kf+''^ w {rf e'^^ 
8 

,2 ' 2p(r) 



0(r) 
N{t) 



v6fce 
2 



In 



{1-Kf 



exp 



[l + K) 
'6ke 



l+K 



Kin' 



iw (r) 



K(t) (r) 



where 



w (t) 
p{t) 



E 



y/iuj^AB - E'^ 



Be- 



: arctan 



2 A'2ujT 



2uj^Ae 



E 



\/4uj-^AB - E^ 



and the constants A, B,pq, C, E,uj are related with the constraint 

— = ~2uj^AB -^-\ -. 

w2 2 2tj2 

Subcase K > 1 In this case the analytic solution of the system of equations (|50ll49p is 

8 

In ( ^i^-Ci U X In + 2p- {r)\ 



N{t) = 



V&ke 
2 



exp 



,(1 + ^) 
^6fce 



l+K 



3w (t) 



X0 (r) 



where 



w (r) 
p(r) 



Ae2 



E 



Be- 



: arctan 



f2uj^Ae^"^* - E 



V4Q^AB~W V V^oj^AB - E^ 

and the constants A, B,po, C, E,uj are related with the constraint 



2C 



^ = -2oj'AB + 



Pi , E^ 



2Q^ 



(69) 
(70) 

(71) 

(72) 
(73) 

(74) 

(75) 
(76) 

(77) 

(78) 
(79) 

(80) 



Our solution reduces to that of [27] if we set -E = 0, that is ^j^dm = and to that of jTl] if we set C = 0, 
that is pb = Mb — 0- Finally if we set C = 0, = we obtain the solution of [22] . 



6 Conclusion 

In order to model all types of matter in the universe we have considered a mixture of three fluids, a perfect fluid 
with a linear EoS, a dust for DM and a scalar field for DE. All fluids are assumed to be minimally interacting 
and the background space is assumed to be the flat FRW spacetime. The available field equations do not 
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suffice to determine the dynamical system. We need two more assumptions, one which will specify the scalar 
field potential and another which will make the system of dynamical equations solvable. In order to define 
these conditions we follow two steps. First we show that the system of the three fluids is equivalent to a two 
dimensional dynamical system moving in under the constraint E ^constant. Then we require that the 
two dimensional system admits an extra Lie symmetry. This requirement fixes the potential to be exponential. 
Requiring further that the solution is invariant under the Lie symmetry we obtain a two parameter family of 
analytic solutions containing a perfect fluid and a (quintessence of phantom) scalar field with an exponential 
potential. 

Requiring next that the Lagrangian defining the two dimensional dynamical system admits an extra Noether 
symmetry, we obtain a new analytic solution which involves all three types of matter, that is, stiff matter, scalar 
field and dust. It is interesting (see Appendix) that in this case the dynamical system reduces to the Ermakov 
Pinney dynamical system which is known to be integrable [501 HI] ; hence we are able to obtain the analytic 
solution. This solution includes previously found solutions in [TTJ [HI [57| as special cases. 

Acknowledgments. We would like to thank Dr Basilakos S for comments and suggestions and the referee 
for useful references. This work has been partially supported from ELKE (grant 1112) of the University of 
Athens. 

APPENDIX 

Analytic solutions for K ^ 1 
Subcase K < 1 

We introduce new coordinates w,p with the relations 

u— ^ V — ^ —. 81) 

1-K ' 1 + K ^ ' 

Then the Lagrangian ([54)) becomes 

2 2 2 

where lo'^ = 1 — K^. This is the Lagrangian for the hyperbolic Ermakov Pinney dynamical system which is 
known to be integrable [20l [21] . The equations of motion are: 

4C 1 

w" + wp'^ -uj^w + ——=0 (82) 

p" + -p'w' = (83) 
w 

The Hamiltonian is: 

i^=l '2_1 2^,2_^^2_2C; 1 ^ (84) 



The solution of the system of equations is 



w (t) = \/yle2"^ - + -86-2^^^ (85) 

PQu; f 2u;^Ae^^- - E \ 

p (r) = , arctan , . (86) 

Vic^^AB - E^ y^/Auj^AB -E-^ J 



Using the Hamiltonian we find the constrain 

„2 ^2 



where A,B,pQ are constants. 
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Subcase K > 1 

In this case we apply the transformation 



l-K ' l + K 

and the Lagrangian (|54p becomes 

2 2 2 

where Co — — 1. This is the Lagrangian for the harmonic Ermakov Pinney dynamical system. 
The solution of the Euler-Lagrange equations is 



w (r) = y Ae^*"^^ ~ ^ + Be-'^i'^'' 

, ipQUj ( 2up- Ae^'-'^'^ - E 

V (t) = , „ ^ arctan - ^ 

V4:Lj^AB - E'^ V VioJ^AB - E^ 

where A, B,po are constants which are related with the constraint equation 

2 2w^ 
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